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Abstract 

The peak forces exerted on soft and rigid samples by a force microscope have been modeled by performing numerical simulations 
of the tip motion in liquid. The forces are obtained by using two contact mechanics models, Hertz and Tatara. We present a com- 
parison between the numerical simulations and three analytical models for a wide variety of probe and operational parameters. In 
general, the forces derived from analytical expressions are not in good quantitative agreement with the simulations when the Young 
modulus and the set-point amplitude are varied. The only exception is the parametrized approximation that matches the results 
given by Hertz contact mechanics for soft materials and small free amplitudes. We also study the elastic deformation of the sample 
as a function of the imaging conditions for materials with a Young modulus between 25 MPa and 2 GPa. High lateral resolution 
images are predicted by using both small free amplitudes (less than 2 nm for soft materials) and high set-point amplitudes. 



Introduction 

The high-resolution imaging of heterogeneous materials, in par- 
ticular soft materials in liquid, by amplitude modulation atomic 
force microscopy (AM-AFM) is an active area of research in 
nanotechnology [1-11]. In AM-AFM, a sharp tip is attached at 
the end of a microcantilever that oscillates at or near its reso- 
nant frequency. When the tip is in close proximity to the 
sample, the amplitude and the phase shift of the oscillation 
change with the strength of the interaction force. The determin- 



ation of the tip-sample interaction force is a major issue in 
dynamic AFM because the force gives access to the materials 
properties of the sample; nonetheless the force is not a direct 
observable. Therefore, several methods have been proposed to 
reconstruct the force in dynamic AFM [12-18]. However, the 
use of force inversion methods has not been generalized in 
AM-AFM because the accuracy of some of the above methods 
is still under study. On the other hand, numerical simulations 
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have been used to determine the maximum repulsive inter- 
action forces, which are referred to as peak forces hereafter [19- 
23]. An analytical scaling law has been deduced to calculate the 
peak forces in air [21]. This method has been applied to deter- 
mine the force on viral capsids in liquid [24]. However, the 
above expressions are often constrained to a specific inter- 
action force model, such as Hertzian mechanics and thus their 
application range is somehow limited. 

Numerical simulations have supported the development of 
AM-AFM by predicting several properties of the tip motion 
[25-27]. Those simulations provide the standards against which 
new experimental or analytical methods should be compared 
[28]. Recently, we have provided a broader numerical insight 
into the interaction forces in AM-AFM [19] by considering 
elastic, viscoelastic, electrostatic double layer and van der 
Waals interactions. 

Here, we perform an extensive computational study of 
AM-AFM to obtain the peak forces of soft (50 MPa) and 
relatively rigid (2 GPa) materials for two different models of 
contact mechanics, namely Hertz [29] and Tatara [30-32]. We 
also provide a comparison between the numerical simulations 
and three analytical expressions [21,33,34]. The dependence of 
the peak force on a wide range of tip-microcantilever prop- 
erties, operational parameters and mechanical properties 
of the sample is analyzed. The Young modulus (E s ) ranges 
from 25 MPa to 2 GPa; the tip radius (R t ) is varied between 5 
and 10 nm; the free amplitude (Aq) goes from 1 to 10 nm and 
the set-point amplitude (A S p) is within the 0.65y4q to 0.95^4q 
range. 

The numerical results are compared to three analytical models, 
the parametrized [21], the average [33] and the linear one [34]. 
The numerical simulations show significant differences from 
the results given by the analytical approximations, although the 
parametrized expression is in good agreement with the Hertzian 
mechanics. The average model follows the trend of the Tatara 
model for the peak forces when varying the set-point amplitude 
for soft samples. For soft materials, the indentation of the tip 
could be higher than A sp . Thus the tip and the sample are in 
permanent contact during the whole oscillation cycle. In fact, 
the ability of exerting small forces and imaging materials in a 
non-invasive manner can be jeopardized because of the effect of 
a static deflection component when A^IAq decreases. We have 
also studied the relationship among peak forces, lateral resolu- 
tion and sample properties for soft (50 MPa) and rigid (2 GPa) 
samples. We deduce a rule to image soft materials with a lateral 
resolution below 3 nm that involves the application of forces in 
the sub- 100 pN regime, the use of cantilevers with force 
constants below 0.1 N/m, free amplitudes below 2 nm and rela- 



tive sharp tips (R t < 5nm). AM-AFM operation at relatively 
high amplitudes can also lead to tip blunting [35,36]. The esti- 
mation of the peak force prior to performing the experiment 
could prevent tip damage. 

Results and Discussion 

Tip motion and contact time for soft and 

relatively rigid materials 

In AM-AFM the equation of motion for the microcantilever-tip 
system is approximated by using the point-mass model [25], 

mz(t) = -kz(t) - z(f) + F ts (d) + F 0 cosotf (i) 

where m is the effective cantilever mass that includes the added 
mass of the fluid, and coq, Q, k and F ts are, respectively, angular 
resonant frequency, quality factor, spring constant and 
tip-sample interaction force. The point-mass model is suitable 
if the contribution of higher modes to the cantilever motion is 
negligible [37]. This could be the case in liquid for small free 
amplitudes, say below 1.5 nm [38]. At higher amplitudes, the 
tip-surface force generates higher harmonics components, 
which could lead to the momentary excitation of higher eigen- 
modes, in particular the second eigenmode [7]. To account for 
those effects we also describe the microcantilever-tip system by 
using an extended Euler-Bernoulli equation [39]. This model 
considers the cantilever as a continuous and uniform rectan- 
gular beam under the action of external forces, 



EI- 



t) + a x 



dw(x, t) 

dt 



+ pbh 



d z w(x,t) 



(2) 



where E is the Young modulus of the cantilever, / the area 
moment of inertia, a\ the internal damping coefficient, p the 
mass density; b, h and L are, respectively, the width, height and 
length of the cantilever; ao is the hydrodynamic damping; w(x,t) 
is the time dependent vertical displacement of the differential 
element of the beam placed at the x position, and F ts tip-sample 
interaction force. 

Equations 1 and 2 are numerically solved by using a fourth- 
order Runge-Kutta algorithm [40]. One should note that the use 
of Equations 1 and 2 in environments of low Q are valid for 
directly excited cantilevers, such as magnetic [41-43] or 
photothermal excitations [44,45]. The tip-sample interaction 
forces are modelled by using two different contact mechanics 
models, Hertz [29] and Tatara [30-32]. The widely used Hertz 
model gives the force as 
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Hertz 



-eff 



3/2 



(3) 



The Tatara contact mechanics has two key differences with 
respect to Hertzian mechanics. First it includes the finite size of 
the sample and second it also considers that the sample defor- 
mation happens symmetrically at both the tip-sample and the 
sample-substrate interfaces. Thus the vertical and lateral 
displacements are part of the contact force computed with this 
model. 



Tatara 



-,3/2 



5 3/2 + 3^ 5 2 + 
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(4) 



where 
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eff 



(5) 



(6) 



(7) 



(8) 



The subindexes "t" and "s" stand, respectively, for tip and 
sample. In the above equations, 5 is the indentation depth, v is 
the Poisson coefficient (v t = 0.3 and v s = 0.4) and E the Young 
modulus with E t = 170 GPa. Each contact mechanics model is 
particularly suited for some types of deformations or strains. 



Hertz contact mechanics is used to deal with small deforma- 
tions at the region of contact between tip and sample. Large 
deformations and finite object sizes are not well described by 
Hertz contact mechanics [32]. As a consequence, the Tatara 
model allows a maximum vertical deformation equal to R Q ff and 
is particularly suited to describe large deformations (with 
respect to the original size) of relatively soft matter, in which a 
vertical force generates both vertical and lateral deformations. 
To apply contact mechanics models in conditions that do not 
meet the model assumptions will lead to unadequate numerical 
estimations. 

We have not found significant differences in the calculation of 
peak forces by using the point-mass model and the continuous 
beam for free amplitudes below 2 nm. For that reason, the data 
for Aq = 1 nm has been obtained with the point-mass 
model while for Aq = 10 nm we have used the extended 
Euler-Bernoulli model. 

Figure 1 shows one period of the tip oscillation and the corres- 
ponding force. The peak force is defined as the maximum force 
point in the dashed line curves. The curves show a purely repul- 
sive interaction, which starts as soon as the mechanical contact 
is established. The tip-sample interface according to Tatara 
(Figure la) or Hertz (Figure lb) is also shown. Both contact 
mechancis models have been applied to describe the response of 
soft (50 MPa) and relatively rigid (2 GPa) surfaces. Figure SI 
of Supporting Information File 1 shows the instantaneous force 
for a variable set-point amplitude. 

Simulated and analytical peak forces values: 
Dependence on the Young modulus and the 
set-point amplitude 

Hu and Raman [21], Kowalewski and Legleiter [34] and 
Rodriguez and Garcia [33] have derived some analytical scaling 
laws to determine the interaction forces in AM-AFM. Hu and 




Figure 1 : Peak forces, tip motion and contact time for two materials, (a) Soft sample (E s = 50 MPa) simulated with the Tatara interaction force (Equa- 
tion 4). (b) Rigid sample (E s = 2 GPa) simulated with the Hertz interaction force (Equation 3). Simulation inputs: k= 0.1 N/m, fo = 25 kHz, Q = 2, 
Aq = 1 nm, fi t = 5 nm, R s = 4 nm and A sp = 0.9 Aq. 
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Raman parametrized the peak force (repulsive) by using a 
nonlinear asymptotic theory [46] and Hertz contact mechanics, 



F =2 
1 par ^ 



1/8 





1/4 r 


71 




T 




V J 





Q 



3 ^effV^t 



1/4 



? 1 \3/8 



Rodriguez and Garcia, by using the virial-dissipation method 
[39,47,48], deduced the following expression for the mean 
value of the force during an oscillation, 



average 



22 



1- 



4T 



-All 



(10) 



some conditions. The parametrized model overestimates the 
peak forces in the case of very soft materials with a maximum 
error of 11%. On the other hand, for stiffer materials the force is 
underestimated with a maximum deviation of 16%. Numerical 
simulations performed with the Tatara model give smaller peak 
force values than those obtained from the Hertz model [19,49]. 
This is because in Tatara contact mechanics the deformation 
happens at both the tip-sample and the sample-substrate inter- 
faces. The linear and average expressions fail to capture the 
trend of the numerical simulations because those expressions 
have been exclusively deduced from the dynamic properties of 
the tip motion and do not consider any influence of the ma- 
terials properties of the sample. Additional comparisons by 
varying the tip radius are presented in Figure S2 of Supporting 
Information File 1 . 



In the absence of long-range attractive forces, the average force 
can provide an estimation of the peak force. 

Kowalewski and Legleiter proposed an extension of the Hooke 
law to determine the force in AM-AFM [34]. In this expression 
the force depends linearly on the amplitude reduction, 

f A^-O (id 

lmear 2Q 

Figure 2 shows a comparison between the analytical expres- 
sions in Equations 9-11 and the numerical results (Hertz and 
Tatara). The comparison covers a range of the Young modulus 
from 25 to 2000 MPa. Hu and Raman parametrized force 
matches the numerical results derived from the Hertz model in 



In Figure 3, the reduction of A sp from 0.95Aq to 0.65Aq 
produces an increase of the peak force. This trend is reported by 
all the approximations and simulations. However, the linear 
approximation give values that are smaller by a factor of 5-100 
compared with the numerical simulations. As a consequence, 
the linear approximation should not be used to estimate the peak 
force in AM-AFM. The average model gives values close to the 
Tatara model for soft materials (25-50 MPa). However, it fails 
to reproduce the data for stiffer surfaces. The average model 
gives the mean value of the forces, attractive and repulsive, 
acting on the tip during an oscillation period. Consequently, 
whenever the forces change significantly with the distance (stiff 
materials) this approximation will fail to give a good estimation 
of the peak force. The parametrized model gives a good numer- 
ical description of the peak forces derived from the Hertz model 
for relatively soft materials. 





Figure 2: Peak force dependence on the Young modulus of the sample for different numerical simulations (Hertz and Tatara) and force analytical 
approximations (parametrized, average and linear), (a) Aq = 1 nm, fi t = 5 nm and (b) Aq = 10 nm, R^ = 5 nm. Other simulation inputs are: k = 0.1 N/m, 
f 0 = 25 kHz, Q = 2, R s = 4 nm, and /A sp = 0.9 A 0 . 
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Figure 3: Dependence of the peak force on the set-point amplitude for different numerical simulations (Hertz and Tatara) and force analytical approxi- 
mations (parametrized, average and linear), (a) Aq = 1 nm and E s = 50 MPa. (b) Aq = 1 nm and E s = 2 GPa. Simulation inputs are: k = 0.1 N/m, 
f 0 = 25 kHz, Q = 2, fi t = 5 nm, and fi s = 4 nm. 



Sample deformation in terms of £ s , Aq, 
and R{ 

The deformation (indentation) exerted by the tip can be consid- 
ered as an indicator of the degree of invasiveness of the tech- 
nique. The dependence of the indentation on the Young 
modulus for Hertz and Tatara models is shown in Figure 4. The 
indentation values are computed and normalized by the free 
amplitude for two values, 1 and 10 nm, and two tip radii, 5 and 
10 nm, respectively. The two-colour curve separates the opera- 
tional parameters in which the deformation is smaller than the 
set-point amplitude from those in which the deformation is 
larger. As expected, the indentation increases by decreasing the 
Young modulus of the sample. Remarkably, for soft materials 
(i.e., those with E s < 100 MPa) the indentation values are close 
to or even larger than the set-point amplitude. This means that 
the tip and the sample are in permanent contact during the 
whole oscillation. This result was observed experimentally by 
Raman et al. [6] while imaging cells. Hertz contact mechanics 
gives smaller indentations than Tatara. In addition, these results 
underline the relevance of the contribution from the static 
deflection, which cannot be neglected in liquid while imaging 
soft materials [19]. However, the above effect decreases, for the 
same ratio A sp /Aq, when increasing the free amplitude as shown 
in Figure 4b. 

Lateral resolution at small peak forces 

Imaging at high-spatial resolution demands a compromise 
between probe, operational parameters and sample properties. 
Figure 5 shows the lateral resolution as given by Hertz and 
Tatara models for two materials, respectively, E s = 50 MPa and 
2 GPa. The lateral resolution is defined as the contact diameter 
between tip and sample. In Figure 5 we visualize the interplay 
between the free oscillation amplitude and set-point amplitude 
with the lateral resolution and peak forces. For a fixed A s JAq 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

E(GPa) 




E s (GPa) 

Figure 4: Normalized indentation as a function of the Young modulus 
of the sample for the Hertz and Tatara models. The indentation was 
normalized to the set-point amplitude: (a) Aq - 1 nm and (b) 
/Aq = 10 nm. In the data above the dashed line indicates the point from 
which the whole oscillation is performed in contact to the material. 
Simulation inputs: k = 0.1 N/m, f 0 = 25 kHz, Q = 2, A sp = 0.65A 0 , 
R s = 4 nm, and two different fi t of 5 and 10 nm, respectively. 

ratio the contact diameter increases with Aq, which reduces the 
lateral resolution. Lowering the A^/Aq ratio down to the range 
between 0.65 and 0.95 also reduces the lateral resolution. In any 
situation the Tatara model gives a better lateral resolution than 
the Hertz model. This result can be traced back to the observa- 
tion that, for the same operational conditions and probe values, 
the Tatara model gives smaller peak forces than the Hertz 
model. The lateral resolution also depends on the elastic 
response of the sample. As a general rule, the stiff er the sample 
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Figure 5: Lateral resolution maps for Hertz and Tatara contact mechanics, (a) E s = 50 MPa. (b) E s = 2 GPa. Filled symbols for >A S p = 0.9/Aq; empty 
symbols for>A S p = 0.7>Ao. Circles for Tatara and triangles for Hertz. Simulation inputs: k = 0.1 N/m, fo = 25 kHz, Q = 2, >A sp = 0.65>4o, Rt = 5 nm, 
R s = 4 nm, and two different A sp of 0.9 and 0.7 Aq, respectively. 



the better the lateral resolution. Sub-nanometric resolution can 
be achieved by using small Aq and maintaining a relatively high 
A S p/Ao ratio for soft materials and rigid materials. It has been 
reported that, in some special situations, also for low A^IAq 
ratios a high resolution can obtained experimentally [28,50]. 
We note that for soft materials (E s = 50 MPa) and in the best 
case scenario (Tatara model) a lateral resolution below 1 nm 
could only be reached by using a free amplitude below 0.5 nm. 
We have separated the plots into regions, for soft materials 
(Figure 5a) a small sub- 100 pN force value is used; while for 
rigid materials (Figure 5b) a sub-1 nN reference value is consid- 
ered. 

Conclusion 

The numerical simulation of the tip motion in amplitude modu- 
lation AFM provides a comprehensive description of the factors 
that control the peak force and the lateral resolution in liquid. 
We have simulated the peak force for two contact mechanics 
models, Tatara and Hertz, and we have calculated three analyt- 
ical approximations, linear, average and parametrized. The 
linear approximation fails to describe qualitatively and quantita- 
tively the peak forces. The average model captures the peak 
force behaviour with the operational parameters but the quanti- 
tative agreement is poor. The parametrized model resembles the 
results given by Hertz for soft materials and small free ampli- 
tudes but its quantitative accuracy decreases by increasing the 
Young modulus. The results show that the discrepancy between 
the analytical and calculated values tends to decrease with 
smaller Young moduli and higher ratio A^IAq. The spatial reso- 
lution depends on the operational parameters, the elastic 
response of the sample, the peak force, and the contact 
mechanics model. The conditions to achieve a high spatial reso- 



lution become more demanding for lower Young moduli of the 
samples. A high spatial resolution in liquid requires the use of 
rather small oscillation amplitudes. Sub-1 nm lateral resolu- 
tions for a soft material of a Young modulus of 50 MPa will 
require the use of a free amplitude of 0.5 nm or less. Lowering 
the free amplitude of the oscillation improves the lateral resolu- 
tion in liquid. The resolution increases in line with the Young 
modulus of the sample, while keeping the operational parame- 
ters constant. The lateral resolution depends on the contact 
mechanics model used to characterize the sample deformation. 
In the Tatara model the sample is finite, consequently the stress 
is relaxed both vertically and laterally, which, for the same 
indentation, provides smaller forces and consequently a better 
resolution that the result given by Hertz model. 

The results presented here provide a good estimation of the 
peak force values experienced by the samples observed with an 
AFM in liquid. However, the simulations have been performed 
without considering hydration layers or viscoelastic effects that 
arise either from the sample or the hydration layer. Those 
effects could modify the peak force values reported here, 
although we do not expect significant changes for the data 
acquired under the conditions for a high spatial resolution (sub- 
5 nm). 

Supporting Information 

Supporting Information File 1 

Additional experimental details. 

[http ://www.beilstein-j ournals.org/bjnano/content/ 

supplementary/2 1 90-4286-4-96-S 1 .pdf] 
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